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ABSTRACT 


The  Auxiliary  Tracking  System  (ATS)  is  being  developed  to 
allow  for  an  accurate  underwater  tracking  <^apability  that  can 
be  deployed  anywhere  in  the  world.  A  simulation  was  created 
to  determine  the  impact  upon  the  relative  accuracy  of  the 
system  due  to  a  normally  distributed  noise.  The  simulation 
was  used  to  look  at  the  impact  of  range  between  targets  and 
the  depth  spacing  of  the  suspended  transducers  on  the  relative 
tracking  accuracy . 


iii 


I««e99lon  For  ^ 

»TIS  GRAil 

DTIC  TA3 

c 

Utif>;>t5riunccd 

a 

1  icotlon— 

By - - — 

AvaliftbUttv  Cedes 

^  LVt/ii  avA/ar 
IpiBt  1  S}>scJal 


4 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . . .  1 

II.  ATS  CONCEPT  .  .............  2 

A.  GENERAL  DESCRIPTION  . 2 

1.  Buoy  System . 3 

2.  Computer  Site .  3 

B.  SYSTEM  OPERATION  . . 4 


III.  SIMULATION  DESCRIPTION  . 

A.  PROGRAM  DESCRIPTION . .  .  . 

1.  Simulation  Inputs  ....... 

2.  Solving  for  Times  . 

3.  /idding  Noise . . 

4.  Siolving  for  Transducer  Locations 

5.  Solving  for  Target  Locations  .  . 

6.  Collecting  Statistics  .  .  .  .  . 

B.  SIMULATIONS  CONDUCTED  ....... 

1.  Noise  Values  ....  . 

2.  Target  Locations  . 

3.  Transducer  Depths  .  .  .  .  .  .  . 

iv 


6 

6 

6 

7 

9 

10 

11 

14 

14 

15 
15 
15 


.  IV.  DATA  ANALYSIS . , . .  16 

V.  CONCLUSIONS  AND  RECOMMENDATIONS  .  18 

A.  CONCLUSIONS  .  . . 18 

B.  RECOMMENDATIONS . 18 

APPENDIX  A  . . •  •  •  20 

APPENDIX  B . 29 

LIST  OF  REFERENCES  . . 37 

BIBLIOGRAPHY  ......  .  38 

INITIAL  DISTRIBUTION  LIST 


39 


t 


I.  IKTRODaCTION 

Within  the  last  several  years,  the  >^avy  has  seen  a 
tightening  of  the  purse  strings  from  congress,  and  as  a 
result,  must  find  more  cost  effective  ';ays  to  perform  the 
veirious  missions  that  are  still  required.  One  area  that  has 
long  been  a  concern  is  the  limited  number  of  underwater 
tracking  ranges.  Alternatives  to  the  development  of  another 
fixed  range  must  be  considered  because  of  the  financial 
impact. 

On  the  west  coast,  there  are  currently  three  areas  where 
underwater  tracking  ranges  are  available.  Washington,  San 
Clemente,  California,  and  Hawaii.  With  the  limited  number  of 
sites  and  the  other  demands  on  range  time  for  research  and 
development,  fleet  ASW  training  opportunities  are  limited. 
Added  to  this  is  the  fact  that  for  many  exercises,  the  units 
must  transit  to  Hawaii  simply  because  that  is  the  only  range 
available.  The  portable  Auxiliary  Tracking  System  (ATS)  was 
conceived  for  these  reasons. 

The  ATS  has  been  designed  as  a  truly  portable  underwater 
tracking  range.  This  thesis  will  look  at  the  theoretical 
tracking  accuracy  of  the  system  currently  under  engineering 
development  and  testing  at  the  Naval  Undersea  Warfare 
Engineering  Station  (NUWES) . 
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II.  ATS  CONCEPT 


A.  GENERAL  DESCRIPTION 

With  the  recent  improvements  made  in  electronics  and  the 
Global  Positioning  Satellite  System  (GPSS) ,  the  dream  of  a 
portable  tracking  range  now  seems  possible  with  current 
technology.  The  engineering  development  model  consists  of 
four  underwater  tracking  buoys  and  a  portable  computer  site. 
The  buoys  are  deployed  either  from  a  ship  or  helicopter  and 
are  recpverable.  For  the  four  buoy  system,  deployment  will  be 
similar  to  that  shown  in  Figure  l.  The  final  system  will 
consist  of  nine  buoys  to  be  deployed  in  various  patterns  based 
upon  the  needs;  of  tlxe  system  being  tested. 


Figure  1.  ATS  Deployment  Configiuration 


1.  Buoy  Systom 

The  buoys  used  in  the  ATS  will  be  free  floating,  self 
contained  buoys.  They  will  have  a  transducer  suspended  below 
at  a  predetermined  depth  (between  50  and  300  meters)  used  both 
for  tracking  the  targets  and  locating  the  other  units  in  the 
tracking  range.  The  acoustic  signal  processor  contained  in 
each  buoy  will  provide  for  rhe  detection  of  the  acoustic  pings 
emitted  by  the  other  sources.  The  GPSS  receiver  in  the  buoy 
will  provide  the  Latitude/Longitude  and  a  very  accurate  time 
standard  for  the  system.  The  on-board  computer  will  control 
the  GPSS  receiver  and  acoustic  signal  processor,  monitor  the 
status  of  the  buoy,  perform  a  system  up/down  self  test,  time 
stamp  and  format  the  position  data  collected  by  the  GPSS 
receiver  and  acoustic  signal  processor  and  send  the  data  to 
the  computer  site  via  an  RF  data  link. 

2 .  Computer  Site 

The  control  station  for  the  ATS  range  will  be  the 
computer  site,  which  will  be  positioned  on  a  ship  located  off 
the  range,  or  on  land  if  the  range  is  set  up  in  coastal 
waters.  The  computer  site  will  be  contained  in  a  portable  van 
allowing  transport  anywhere  in  the  world  either  by  ship  or  by 
ceurgo  aircraft.  The  only  support  needed  will  be  AC  power,  as 
the  computer  site  will  be  self  contained.  The  computer  site 
will  contain  the  RF  data  link  receivers,  RF  voice 
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communications  equipment,  a  GPSS  receiver,  t ■  range  tracking 
computers  and  peripherals,  and  other  eqo.ipm'jnt  as  needed. 

B.  SYSTEM  OPERATION 

The  concept  embodied  by  the  ATS  is  a  "sing-around”  self 
calibrating  tracking  range.  This  means  that  the  buoys  will 
locate  each  other.  Each  buoy  has  an  acoustic  pinger  on  the 
floating  buoy  and  a  transducer  suspended  below.  In  a  two 
minute  cycle,  each  buoy  pinger  will  transmit  an  acoustic 
signal  followed  by  each  suspended  transducer  also  transmitting 
an  acoustic  signal.  The  acoustic  processor  within  each  buoy 
will  control  the  pinging  cycle;  a  coded  pulse  will  be  used  to 
identify  the  signal.  The  transit  time  of  the  signal  will  be 
known  as  each  buoy  will  transmit  a  pulse  at  a  specific  time, 
and  the  reception  time  of  the  signal  wili  be  sent  to  the 
computer  site  via  the  RF  link  from  the  buoy  that  received  the 
signal.  From  this,  the  transit  time  of  the  signal  is  known, 
and  the  range  between  the  source  and  receiver  is  calculated  by 
using  the  Sound  Velocity  Profile  of  the  water  mass.  The 
signal  from  thf^  buoy  pinger  to  the  transducer  suspended  below 
is  used  to  calculate  the  depth  of  each  transducer.  The 
transit  times  from  transducer  to  transducer  are  then  used  to 
calculate  the  slant  ranges  when  the  depth  of  each  transducer 
has  been  determined.  From  this  data,  a  range  grid  is  set  up 
as  shown  in  Figure  2.  The  location  of  each  transducer  and  the 
range  grid  will  be  updated  every  two  minutes  by  this  method. 
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Flgura  i,  ATS  Rang*  Grid 

Vitti  a  relative  grid  established,  the  system  now  determines 
the  positions  of  the  targets  on  the  range.  The  ATS  has  three 
distinct  frequency  bands  to  utilise,  by  incorporating  a  Shift 
Frequency,  Shift  Keying  (SFSK)  processor  and  transmitter.  The 
high  frequency  will  be  used  for  the  sing-around  calibration 
already  described,  and  the  mid  and  low  frequencies  will  be 
used  to  track  the  torpedo  and  the  target.  The  SFSK  processor 
and  transmitter  is  similar  to  the  system  used  on  the  Quinalt 
Shallow  Water  range  at  NUWFS. 
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III.  SIMULATION  DESCRIPTION 

A.  PROGRAM  DESCRIPTION 

A  simulation  program  was  written  to  determine  the  effects 
of  a  normal  error  on  the  transit  times  received  at  the 
computer  site  to  investigate  the  accuracy  of  the  ATS.  The 
effects  of  shadow  zones  were  not  considered  in  the  simulation 
due  to  the  relatively  short  distances  between  any  source  and 
receiver  pair.  The  simulation  program  has  three  main 
sections;  solving  for  transit  times,  solving  for  positions, 
and  collecting  statistics. 

1.  Simulation  Inputs 

The  ATs  simulation  program  requires  three  sets  of 
inputs.  The  first  set  of  inputs  is  the  sound  speed  profile 
for  the  water  mass.  The  program  can  accept  up  to  60  data 
points  for  the  profile.  The  second  set  of  inputs  is  the 
actual  positions  for  the  four  buoys,  the  four  suspended 
transducers  and  the  two  targets.  The  inputs  are  cartesian  co> 
ordinate  with  buoy  one  being  at  the  origin  and  buoy  three 
establishing  the  X*Z  plane.  The  third  sat  of  inputs  is  for 
Che  number  of  replications  desired  and  the  standard  deviation 
of  the  noise  for  the  systems .  The  standard  deviation  is  a 
function  of  the  distance  between  the  source  and  receiver,  and 
within  the  program  is  modeled  by  changing  the  travel  times  in 


6 


each  of  three  groups ;  buoy  to  transducer ,  transducer  to 
transducer,  and  target  to  transducer. 

2.  Solving  for  Times 

Once  the  inputs  have  been  entered,  the  program  solves 
for  the  actual  transit  time  between  each  source  and  receiver 
pair.  This  is  done  by  taking  the  slant  range  between  each 
source  and  receiver  (Equation  3.1)  and  dividing  by  the 

Ranges,! (X,-X^) ^  (3*1) 

effective  sound  speed  between  source  and  receiver.  The 
effective  sound  speed  is  calculated  jy  solving  for  the  travel 
time  for  a  ve^rtically  propagating  wave  in  each  layer.  The 
travel  times  are  summed  between  the  source  and  receiver  depths 
and  the  difference  in  depths  is  divided  by  the  travel  time  to 
yield  the  effective  sound  speed  ^tween  the  source  and 
receiver.  The  derivation  of  the  equation  is  shown  below  and 
starts  with  the  definition  of  the  travel  time  (Equation  3.2). 

The  ray  path  within  a  constant  gradient  is  defined  as  the  arc 
of  a  circle  of  radius  equal  to  a  (Equation  3.3  and  3.4) 

[Ref  l:p.  402]. 

ds*oeid  (3.3) 
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o=-.£l£L 

g  COS0 


(3.4) 


Substitution  of  Equation  3.4  into  Equation  3.2  yields  Equation 
3.5  and  3.6.  The  solution  of  the  integral  >  is  shown  in 
Equation  3.7 


CO 


COS0 


t*—  f  sec  (6)  d9 


t=-i  [ln|sec(8)  4-tan(8)  |1 


(3.5) 


(3.6) 


(3.7) 


Separation  of  the  terms  yield  Equation  3.8.  Using  snell's  law 
[Ref  1;  p.  401],  as  shown  in  Equation  3.9,  and  making  the 
assumption  that  the  sine  terms  cancel  for  a  vertically 
propagating  wave.  Equation.  3.8  can  be  reduced  to  Equation 
3.10. 


i,„|i>8ine  cose. 
g  cos8  i+sin8_' 


cose^,  C(Zg) 

cos8  c(z) 


(3.8) 


(3.9) 


g  (c(Zji) 


(3.10) 
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By  summing  the  travel  times  within  each  layer,  the 
propagation  time  is  found.  If  within  a  layer,  the  gradient  is 
equal  to  zero,  the  program  takes  the  vertical  distance  within 
the  layer  and  divides  by  the  sound  speed  within  the  layer  to 
determine  the  time  t.  The  effective  sound  speed  is  now  found 
(Equation  3.11). 


effspd= 


(3.11) 


Knowing  the  slant  range  and  the  effective  sound  speed  between 
each  source  and  receiver  pair,  the  transit  times  for  all  pairs 
are  determined  (Equation  3.12). 


effspd 


(3.12) 


3.  Adding  voise 

When  all  of  the  transit  times  have  been  calculated  and 
stored,  the  program  begins  the  iteration  cycle.  The  firdt 
step  is  the  addition  of  the  noise  value  to  the  transit  times. 
This  is  done  by  taking  a  standard  nbimal  random  variable  and 
multiplying  it  by  the  transit  time  and  the  percentage  value 
for  the  standard  deviation.  This  results  in  a  normal  random 
variable  with  mean  equal  to  zero  and  standard  deviation  equal 
to  the  given  percentage  of  the  transit  time.  This  method  was 
chosen  because  the  transducer  to  transducer  ranges  and  times 


are  typically  5  times  greater  than  the  pinger  to  suspended 
transducer  times,  and  the  error  should  be  a  function  of  range. 

4.  Solving  for  Transducer  Locations 

After  the  noise  has  been  added  to  each  of  the  transit 
times,  the  location  of  each  transducer  is  calculated.  The 
depth  of  each  transducer  is  found  first  using  the  transit  time 
from  the  buoy  pinger  to  the  transducer  suspended  below  (i.e. 
buoy  pinger  one  to  transducer  one) .  From  this  transit  time 
and  using  the  length  of  the  cable  as  a  starting  point,  the 
depth  is  determined  through  an  iterative  process  by 
calculating  the  effective  sound  velocity  between  the  pinger 
and  the  estimated  depth  of  the  transducer.  The  solution  is 
then  found  by  multiplying  the  transit  time  by  the  effective 
soxind  velocity.  The  resulting  depth  is  then  compared  with  the 
estimated  depth,  and  when  the  difference  between  the  two 
values  is  less  than  one  foot,  the  program  continues,  otherwise 
the  resulting  depth  is  used  for  the  estimated  depth  and  the 
calculation  is  repeated.  Once  the  four  transducer  depths  are 
established,  the  range  coordinate  system  is  established  using 
transducer  one  as  the  origin  of  the  X-Y  grid  and  transducer 
three  as  the  intersection  with  the  horizontal  axis  (X-axis) . 
There  are  two  transit  times  for  each  range  calculation  made 
(transit  time  from  transduceir  one  to  t-'iree  and  from  transducer 
three  to  one) ,  and  the  times  are  averaged.  The  slant  range 
(SR)  is  then  found  by  calculating  the  effective  scund  velocity 
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between  the  two  depths  and  multiplying  by  the  average  transit 
time  between  the  two  transducers.  The  horizontal  distance  (H) 
is  calculated  using  Equation  3.13.  The  position  of  transducer 

(3.13) 

two  is  found  by  using  the  horizontal  ranges  from  transducers 
one  and  three  and  using  the  GPSS  position  of  buoy  two  to 
determine  which  of  two  unique  solutions  is  correct  (the 
intersection  of  two  circles) .  Transducer  four  is  then  located 
the  same  way  using  the  positions  and  ranges  from  transducers 
two  and  three. 

S.  Solving  for  Target  Locations 

With  the  range  coordinaite  system  established  and  the 
location  of  all  transducers  found,  we  are  positioned  to 
calculate  the  location  of  the  two  targets.  It  is  assumed  that 
all  four  transducers  receive  the  signal  from  each  target.  The 
location  of  the  target  is  then  calculated  through  a  least 
squares  method  similar  to  one  used  by  Larry  Anderson  at  NUWES. 
The  detailed  procedures  for  the  algorithm  are  presented  below. 

Given  the  location  of  each  transducer  as  (Xj,  Yj,  ZJ  and 
the  solution  range  between  transducer  i  and  the  target  (R.^)^ 
from  the  transit  time,  the  location  of  the  Target  (Xf,  \’j,  Z^) 
is  calculated.  is  the  measured  slant  range  between  target 
and  transducer  i,  and  R|  is  the  calculated  slant  range  between 
target  and  transducer  1  as  given  by  Equation  3.14.  The 
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(Xr-Xi)  2+  ( Yj.-Yi)  (Zy-Zi)  ‘ 


(3.14) 


location  of  the  target  is  then  calculated  with  respect  to 
minimizing  the  Sum  of  Squares  of  errors  (SSE) ,  where  the  SSE 
is  defined  by  Equation  3.15.  The  calculated  slant  ranges  (Rj) 


4 

in 


(3.15) 


are  linearized  around  the  trial  value  (Xt6»  ^To»  /  shown  in 
Equation  3.16,  where  means  evaluated  at  point  '’o'*.  The 


dRi 


(X  -X  ) 


dRi 


(Z^-ZrJ  (3.16) 

3*  A 


partial  derivatives  are  listed  in  Equations  3 . 17  through  3.19. 


dXf 


X.-Xi 


R 


i. 


(3.17) 


dRi 

ay. 


o' 


(3.18) 


dR^ 

az. 


z^ 


(3.19) 


The  linearized  model  is  substituted  into  Equation  3.15  and  the 
SSE  is  minimized  with  respect  to  (X,,  Y„  Z,) .  The  resulting 
normal  equations  are  put  in  matrix  form  as  shown  in  Equation 
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V 


3.20.  The  individual  matrices  are  shown  in  Equations  3.21 
through  3.23.  Equation  3.20  is  solved  for  the  correction 
matrix  (Equation  3.22)  which  when  added  to  (X^o#  Xto*  Zt^) 


AX==B 


(3.20) 


A 


2 

dRi 

_  8J?j 

^R^ 

o 

3iv 

IK 

o  ~ 

S-sz; 

o  * 

"Wrl 

O 

2 

—  SR, 

dR^l 

^Ix, 

dy, 

0  Tl 

o  ^ 

o' 

S-Sz; 

■diV 

o  ^ 

dR^ 

azi 

o  ^ 

r-  3R, 
^  dYr 

O  ^ 

dRi 

'SK 

o  * 

Y-  * 

^  dZr 

o  ^  o 

(3.21) 


X 


Xf-Xf^ 

Yr-Yr, 


(3.22) 


o 

B  • 

o 

o 

(3.23) 


yields  the  improved  target  position  (Xj,  Yj,  Zj) .  If  the  sum 
of  the  squares  of  the  corrections  is  less  than  some  value 
epsilon,  Equation  3.24,  the  iteration  is  stopped.  Otherwise, 


(Xr-Xr)  ( Yj-YrJ  (Zy-Z,,)  *<e 


(3.24) 
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the  target  position  (Xj,  Z,-)  is  used  for  the  next  iteration 
in  place  of  (Xto#  Y^o,  Zto)  and  the  procedure  is  repeated  until 
50  iterations  have  been  completed. 

6.  Collecting  Statistics 

When  the  positions  of  the  targets  have  been  solved 
for,  the  errors  are  calculated.  The  distance  between  the 
actual  and  solution  location  of  each  target  (Equation  3.25), 
and  the  distance  in  the  relative  position  of  the  targets  with 
respect  to  each  other  (Equation  ,3.26)  are  calculated. 

(Zol-Zal) " 

)^*...  (3.26) 

X,j  is  the  X-value  for  the  true  position  of  the  ith  target,  X,i 
is  the  X-value  for  the  solved  position  of  the  ith  target  ,  TE, 
is  the  True  error  for  the  ith  target,  and  RE  is  the  relative 
error  between  the  two  targets.  The  values  for  TEj  and  RE  are 
summed  over  the  iteration  loop  and  the  programs  final  output 
is  the  mean  and  standard  deviation  of  the  true  position  error 
for  each  target  and  the  relative  position  error  between  the 
two  targets. 

B.  BIMULATIONS  CONDUCTEB 

Three  factors  are  considered  as  to  the  influence  of  the 
relative  accuracy:  Magnitude  of  the  noise  in  the 
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measursments ,  the  proximity  of  target  one  to  target  two,  and 
whether  the  transducers  are  at  equal  or  staggered  depths. 

1.  Noise  Values 

Noise  values  are  entered  into  the  simulation  by  the 
standard  deviation  being  a  percentage  of  the  true  value. 
Values  examined  range  between  0.1  and  1.0  percent. 

2.  Target  Locations 

The  simulations  were  conducted  for  targets  at  close 
range  (100  ft  separation)  and  long  range  (3000  ft  separation) 
to  investigate  the  effects  of  target  separation  on  the 
relative  accuracy. 

3.  Transducer  Depths 

Some  discussion  has  taken  place  as  to  the  effects  of 
staggering  the  transducer  depths  or  having  them  all  at  an 
equal  depth  with  the  thought  being  would  accuracy  be  gained  by 
having  the  suspended  transducers  staggered  at  varying  depths. 
Would  this  allow  for  better  target  resolution  given  the  four 
slant  ranges  from  the  target  to  the  transducers?  Simulations 
were  also  conducted  to  draw  a  conclusion  on  the  selection  of 
transducers  depths;  equal  depths  or  staggered  depths  with 
intervals  of  100  feet. 
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IV.  DATA  AMALYSZS 


Each  of  the  four  following  cases  were  run  using  the 
simulation  program  with  the  noise  values  ranging  from  0  to  1 
percent.  The  performance  goal  for  the  ATS  is  20  yards 
relative  accuracy.  The  specifics  of  the  four  cases  are  as 
follows: 

1.  Case  One:  Transducers  at  equal  depth  (1000  feet)  and 
targets  at  close  range  (100  feet). 

2.  Case  Two?  Transducers  at  equal  depth  and  targets  at 
long  range  (3000  feet) . 

3.  Case  Three:  Transducers  at  staggered  depths  (100  foot 
intervals)  and  targets  at  close  range. 

4.  Case  Four:  Transducers  at  staggered  depths  and  targets 
at  long  range. 

The  results  for  the  four  cases  are  shown  in  Figure  3.  The 
limit  for  the  noise  varies  from  0.3  percent  for  cases  two  and 
four  (targets  at  long  range)  to  0.4  percent  for  cases  one  and 
three  (targets  at  close  range) . 


Case  1  Case  2  Case  3  ••«•••  Case  4 


Figure  3.  Relative  Error  Plot 
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V.  CONCLUSIOKS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  results  from  the  simulation  show  that  the  system  can 
not  handle  a  lot  of  noise  as  currently  designed.  The 
staggering  of  the  transducer  depths  did  not  appear  to  have  an 
effect  on  the  accuracy  of  the  system.  All  four  cases  were 
relatively  close  for  noise  value  where  the  relative  accuracy 
exceeded  20  yards. 


B.  RECOMMENDATIONS 

The  accuracy  assessment  for  the  ATS  is  a  very  complex 
problem.  This  study  has  only  focused  on  a  few  -areas  and  a 
complete  assessment  of  the  system  will  require  further  study. 
From  the  initial  results  on  the  engineering  development  model, 
the  outlook  is  very  promising  for  the  Navy  to  receive 


substantial  benefits  from  its  continued  funding. 


the  data  from  the  prototype  buoys  should  reveal  the  magnitude 
of  the  noise,  and  the  expected  accuracy  can  then  be 
determined.  An  additional  area  that  should  also  tie  considered 
is  to  use  telemetry  data  from  the  targets  in  the  tracking 
problem.  This  could  be  utilized  to  provide  the 
target  and  the  velocity.  By  incorporating  a  Kklman  Filter 
into  the  algorithm,  it  may  be  possible  to  reduca  the  impact 


IS 


Analysis  of 


from  a  bad  ping  cycle.  This  would  also  necessitate  looking  at 
the  accuracy  of  the  depth  sensors  currently  available. 

The  Auxiliary  Tracking  System  should  continue  to  be 
developed.  Potential  benefits  for  the  Navy  warrant  the 
investment  now,  as  current  ranges  are  taxed  to  keep  up  with 
the  demands  placed  upon  them  by  the  many  users.  As  new 
weapons  systems  are  developed,  the  ATS  will  allow  more 
realistic  operational  testing  to  be  conducted  by  allowing  the 
weapons  systems  to  be  evaluated  in  various  ocean  types. 
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APPENDIX  A 


CC  ATSSIM.FOR 

CC  SIMULATION  PROGRAM  TO  EVALUATE  AUXILIARY  TRACKING  SYSTEM 

CC  ABSOLUTE  AND  RELATIVE  ACCURACY  ASSESSMENTS  BASED  UPON  A 

CC  RANDOM  NOISE  WIV\  A  GIVEN  MEAN  AND  VARIANCE. 

CC 

REAL*8  SEED 

REAL*4  XBUOY,  YBUOY,ZBUOY,XPHONE,YPHONE,ZPHONE,XTGT,YTGT,ZTGT, 

*  MEANl,  VARl,  MEAN2,  VAR2,  MEANR,  VARR,  SRANGE,  RANGE,  NOISE, 

+  TTIMEl,  HTIMEl,  BTIMEl,  TTIME,  HTIME,  BTIME,  A,  B 
INTEGER* 4  BUOYS,  MAXSSV,  NCOUNT 

DIMENSION  XBUOY ( 4 ) , YBUOY ( 4 ) , ZBUOY ( 4 ) , XPHONE { 4 ) , YPHONE  ( 4 ) , 

+  ZPHONE(4),  SRANGE(6),BTIME1(4,4),HTIME1(4,4), 

+  PHONEX(4),  PHONEY(4),PHONEZ(4),  DEPTH ( 40 ) , SSPVSL { 40 ) , 

+  XTGTr2) ,YTGT(2) ,ZTGT(2) ,TGTX(2),TCTY(2),TTIME1(2,4) , 

+  TGT2(2),BTIME(4,4),HTIME(4,4),TTIME{2,4), 

+  TERROR(2),TESQ(2),TESOM(2),TRANGE(2,4) 

MAXSSV  >  40 
BUOYS  «  4 
ZCABLE  »  1000.0 

CALL  INPUT  (NCOUNT,  BTN,  STN,  VTN,  SEED) 

CALL  PROFILE  (MAXSSV, NUMSSV,DEPTH,SSPVEL) 

CALL  POSITS  ( XBUOY, YBUOY, ZBUOY, XPHONE, YPHONE, ZPHONE, 

+  XTGT,YTGT,ZTGT)  , 

CC  CALCULATE  TRANSIT  TIMES  FOR  BUOYS 
DO  20  I  -  1,  4 
DO  10  J  «  1,  4 

RANGE  »  SQRT( (XBUOY(l)-XPHONE(J) )**2+(YBUOY(I)-YPHONE(J) )**2  + 
+  (ZBDOY(I)  -  ZPHONE ( J) )**2) 

CALL  HMVELOC ( ZBUOY ( I ) , ZPHONE ( J) , NUMSSV, DEPTH , SSPVEL, HHVEL ) 
BTIME ( I, J)  «  RANGE/HMVEL 
10  CONTINUE 
20  CONTINUE 

DO  so  I  »  1,  4 
DO  40  J  •  1,  4 

IF  (1  .EQ.  J)  GOTO  30 
RANGE  •  SQRT((XPH0NE(I)*-XPH0NE(J))**2+ 

+  (YPHONE ( .• ) -YPHONE (J) ) **2+( ZPHONE ( I ) -ZPHONE (J) ) **2 ) 

CALL  HMVELOC(ZPHONE(I) ,  Z?HONE(J),  NUMSSV,0EPTH, SSPVSL, HHVEL) 
HTIME(I,J)  «  RANGE/HMVEL 
30  CONTINUE 
40  CONTINUE 
SC  CONTINUE 

DO  70  I  •  1,  2 
DO  60  J  -  1,  4 

RANGE  ■  SQRT( (XTGT(I)-XPHONE(J) )**2+(YTCT(I)-YPHONE(J) )**2+ 

+  (ZTGT(I)  -  ZPHONE (J) )**2) 

CALL  HMVELOC ( ZTGT ( I ) , ZPHONE ( J ) , NUMSSV , DEPTH , SSPVEL , HMVEL ) 

TTIME ( I, J)  •  RANGE/HMVEL 
60  CONTINUE 
70  CONTINUE 
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DO  900  KOUNT»l,  NCOUNT 

79  CONTINUE 
CC  ADD  NOISE 

DO  71  I  «  1,  2 
J  «  I  +  2 

CALL  URAND(A,B,SEED) 

BTIME1(I,I)  «  ETIME(I,I)  +  (A*BTN*BT1ME{I,I) ) 

BTIME1(J,J)  »  BTIME(J,J)  +  (B*BTN*BTIME( J, J) ) 

71  CONTINUE 

DO  73  I  -  1,  4 
DO  7r  J  -  1,  2 
K  -  J  +  2 

CALL  UBANO  (A, B, SEED) 

HTIME1{I,J)  -  HTIME(I,J)  (A*STN*HTIME(I,  J) ) 

HTIME1(1,K)  »  HTIME(I,K)  +  (B*STN*HTIME(I,K) ) 

72  CONTINUE 

73  CONTINUE 

DO  75  J  »  1,  4 
CALL  URAND(A,B,SEED) 

TTIME1(1,J)  -  TTIME(1,J)  +  (A*TTN*TTIMB(1, J) ) 

TTIME1(2,J)  «  TTIME(2,J)  +  (B*BTN*TTIME(2, J) ) 

75  CONTINUE 

CC  CALC  BUOY  POSITS 
DO  80  I  >  1,  4 
PHONEX(I}  -  XBUOY(I) 

PHONEY(I)  -  YBDOY(I) 

80  CONTINUE 

CALL  HYDEPTH  (ZCABLE,  ZBUOY,  NUMSSV,  DEPTH, 

♦  SSPVEL,  HHVEL,  BTIMEl,  PHONEZ) 

CALL  RPHONE  (SRANGE,  HTIMEl,  PHONEZ,  NUMSSV,  DEPTH,  SSPVEL,  HMVEL) 
CALL  TWOD  (SRANGE,  PHONEX,  PHCSEY,  PHONEZ,  YBUOY) 

CC  CALC  TARGET  POSITS 
DO  100  J  -  1,  2 
TCTZ(J)  -  400 
TCTX(J)  -  2000 
TGTY(J)  »  2000 

CALL  TGTRANGE  (TTIMEl,TGTZ,Z?HONE, NUMSSV, DEPTH, SSPVEL, HMVEL, 

•f  TRANGE.J) 

CALL  THREED  (J,TGTX,TCTY,TGTZ, PHONEX, PHONEY, PHONE Z,TRANGE) 

100  CONTINUE 
CC  CALC  TRUE  ERRORS 
ERFIX  -  0.0 
DO  150  J  -  1,2 

TERROR(J)  •  SQRT((XTCT(J)-TCTX(J))**2+(YTGT(J)-TGTY(J))**2+ 

+  {ZTGT(J)-TGTZ(J))**2) 

IF  (TERROR(J)  .GT.  1000.0)  ERFIX  >  1.0 
150  CONTINUE 

IF  (ERFIX  .EQ.  1)  THEN 
GOTO  79 
ELSE 

DO  200  J  -1,2 

TESQ(J)  -  TESQ(J)  +  TERROR(J)**2 
TESUM(J)  -  TESUM(J)  ♦  TERROR(J) 

200  CONTINUE 

CC  CALC  RELATIVE  ERRORS 

RERROR  -  SQRT((XTCT(1)“TGTX(1)-XTCT(2)+TGTX(2))»*2+ 

+  (YTGT(1)-TCTY(1)-YTCT(2)+TGTY(2) )**2+ 

+  (ZTCT(1)-TGTZ(1)-ZTCT(2)+TCTZ(2) )**2) 

RBSQ  -  RESQ  *  RERROR**2 
RESUM  -  RESUM  RERROR 
ENDIF 
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900  CONTINUE 

CC  CALC  STATS  ON  TRUE  ERRORS 
CC  TGTl  RESULTS 

MEANl  »  TESUM(l)/NCOUNT 

VARI  ■  SQRT(ABS(NCOUNT*TESQ(l)-TESUM(l)**2) )/(NCOUNT*{NCOUNT-rl) ) 
CC  TGT2  RESULTS 

MEAN2  »  TESUM(2)/NCOUNT 

VAR2  -  SQRT(ABS(NCOUNT*TESQ(2)-TESUM(2)**2) )/(NCOONT*(NCOUNT-l) ) 
CC  CALC  STATS  ON  RELATIVE  ERRORS 
CC  RELATIVE  RESULTS 

MEANR  >  RESCM/NCOUNT 

VARR  ■  SQRT(ABS(NCOUNT*RESQ-RESUM**2) )/(NCOUNT*(NCOUNT-l) ) 

CALL  OUTPUT  (MEANl,  VARI,  MEAN2,  VAR2,  MEANR,  VARR, 

+  BTN,STN,TTN,NCOUNT) 

END 


INPUT  SUBROUTINE  TO  ENTER  IN  SIMULATION  VALUES 


SUBROUTINE  INPUT  (NCOUNT,  BTN,  STN,  TTN,  SEED) 
REAL* 8  SEED 

OPEN  (UNIT  -  1,  PILE  -  'INPUT.DAT') 

READ  (1,*)  NCOUNT,  BTN,  STN,  TTN,  SEED 
CLOSE  (UNIT  «  1,  STATUS  -  'KEEP') 

RETURN 

END 


*  PROFILE  SUBROUTINE  TO  ENTER  SOUND  VELOCITY  PROFILE 


SUBROUTINE  PROFILE  (MAXSSV,  NUMSSV,  DEPTH,  SSPVEL) 

REAL*4  DEPTH,  SSPVEL 
INTECER*4  MAXSSV,  NUMSSV 
DIMENSION  DEPTH(40),  SSPVEL(40) 

OPEN  (UNIT  «  1,  FILE  •  'SSPROF.DAT’) 

READ  (1,*)  NUMSSV 
IF  (NUMSSV  .CT.  MAXSSV)  THEN 

WRITE (*,*)  '  MAXIMUM  NUMBER  OP  SSP  DATA  POINTS  EXCEEDED' 
ENDIF 

DO  10  I  -  1,  NUMSSV 

READ  (1,*)  DEPTH(I),  SSPVEL(I) 

10  CONTINUE 

CLOSE  (UNIT  -  1,  STATUS- 'KEEP') 

RETURN 

END 


POSITS  SUBROUTINE  TO  ENTER  BUOY,  TRANSDUCER  AND  TARGET  • 


SUBROUTINE  POSITS  ( XBUOY , YBUOY , ZBUOY , XPHONE , YPHONE , ZPHONB , XTGT , 
♦  YTCT,ZTCT) 

REAL*4  XBUOY , YBUOY , ZBUOY , XPHONE , YPHONE , ZPHONB 
DIMENSION  XBUOY(4),YBUOY(4),ZBUOY(4),XPHONE(4),YPHONE(4), 
ZPHONE(4),XTCT(2),YTCT(2),ZTCT(2) 

OPEN  (UNIT  -  2,  FILE  -  'BUOYL0C.DAT') 

DO  10  I  -  1,  4 

READ(2,*)  XBUOY(I),YBUOY(l),ZBUOY(I),XPHONE(I),YPHONE(I), 
ZPHONE(l) 

10  CONTINOT 

OPEN  (UNIT  -  1,  FILE  -  'TAPCET.DAT') 

DO  20  J  -  1,  2 


READ(1,*)  XTGT(J),  YTCT(J),  2TGT(J) 
20  CONTINUE 
RETURN 
END 


*  HMVELOC  SUBROUTINE  TO  CALCULATE  MEAN  SOUND  VELOCITY 


SUBROOTINB  HMVELOC  (ZPINC,  ZRCVR,  NUMSSV,  DEPTH,  SSPVEL,HMVEL) 
REAL*4  ZPINC,  ZRCVR,  DEPTH,  SSPVEL,  HMVEL,  CRADNT,  HMVl,  HMV2,T 
INTEGBR*4  I,  J,  K,  INIT,  PINAL,  NUMSSV 
DIMENSION  DEPTH ( NUMSSV ) ,  SSPVEL(NUMSSV) 

C  ORDER  THE  DEPTH  OF  THE  RECEIVER  AND  HYDROPHONE  SO  THAT  THE 
C  RECEIVER  IS  THE  DEEPER  OP  THE  TWO  DEPTHS. 

T  «  0.0 
INTE  -  0 

IP  (ZPINC  .CT.  ZRCVR)  tHEM 
TEMP  «  ZPINC 
ZPINC  -  ZRCVR 
ZRCVR  -  TEMP 
INTE  -  1 
KNOIF 

C  INTERPOLATE  THE  SOUND  SPEED  AT  THE  PINCER  DEPTH 
DO  1  I  «  2,  NUMSSV 

IP  (ZPINC  .CT.  OEPTH(I))  GOTO  1 

CRADNT  •  (SSPVEL(I)  -  SSPVBL( I-l) ) / (DEPTH( 1 )  -  DEPTH(I-l)) 
HMVI  ■  SSPVEL(I-l)  ♦  CRADNT* ( ZPINC-DEPTH ( I-I ) ) 

INIT«  I 
CO  TO  2 

1  CONTINUE 

C  INTERPOLATE  THE  SOUND  SPEED  AT  THE  HYDROPHONE  DEPTH 

2  DO  3  J  -  2,  NUMSSV 

IP  (ZRCVR  .CT.  DBPTM(J))  GOTO  3 

CRADNT  •  (SSPVEL(J)  *  SSPVEL(J>1) )/(DEPTH(J)  -  DEPTH(J-l)) 
HMV2  «  SSPVEL(J>1)  ♦  CRADNT* ( ZRCVR-DEPTH(J-1 ) ) 

PINAL  >  J-1 
CO  TO  4 

3  CONTINUE 

C  CALCULATE  THE  HARMONIC  SOUND  SPEED  BETWEEN  THE  PINCER  AND  THE 
C  HYDROPHONE 

4  IP  (PINAL  .LT.  INIT)  THEN 
fc  ' 

C  SOURCE  AND  RECEIVER  IN  SAMS  GRADIENT 

C  TAU  AVERAGE  OP  TWO  SOUND  SPEEDS  OP  SOURCE  AND  RECEIVER 
C 

CRADNT  •  (KMV2  -  HMVI)  /  (ZRCVR  -  ZPINC) 

IP  (CRADNT  .NS.  0.0)  T  >  LOG (HMV2 /HMVI) /CRADNT 
IP  (CRADNT  .EQ.  0.0)  T  •  (ZRCVR  -  ZPINO)/HMV2 
HMVEL  -  T/(ZRCVR-ZPINC) 

SLSSIP  (PINAL  .SQ.  INIT)  THEN 
C 

C  SOURCE  AND  RECEIVER  SPSRATED  BY  ONE  DEPTH  GRADIENT 
C  ' 

CRADNT  ■  (SSPVSL(INIT)  -  HMVI)/(DZPTH( INIT)  -  ZPINC) 

IP  (CRADNT  .NE.  0.0)  T  >  LOO(SSPVEL(INIT)/HMVI)/CRAONT 
IP  (CRADNT  .SQ.  0.0)  T  *  (0EPTH( INIT)oZPINC) /SSPVSL( INIT) 
CRADNT  -  (KNV2-SSPVSL(PINAL))/( ZRCVR  -  DEPTH(PINAL) ) 

IP  (CRADNT  .NS.  0.0)  T  -  T  *  LOG (HMV3 /SSPVEL ( PINAL) ) /CRADNT 
IP  (CRADNT  .EQ.  0.0)  T  •  T  ♦  ( ZRCVR-OEPTH( PINAL) )/HMV2 
HMVEL  •  ABS( (ZRCVR-ZPINC)/T) 

ELSE 
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SOURCE  AND  RECEIVER  SEPERATED  BY  MORE  THAN  ONE  DEPTH 

GRADNT  *  (SSPVEL{INIT)  -  HMVl )/{ DEPTH ( INIT)  -  ZPING) 

IF  (GRADNT  .NE.  0.0»  T  -  LOG {CSPVEL(  INIT) /HMVl) /GRADNT 
IF  (GRADNT  .EQ.  0.0)  T  -  (DEPTH( INIT) -ZPING) /SSPVEL( INIT) 

DO  5  K-  INIT  +  1,  FINAL 

GRADNT  -  (SSPVEL(K)-SSPVEL(K-l)  )/(DEPTH(K)-DEPTH(K-l) ) 

IF  (GRADNT  .NE.  0-0)  T  »  T+LOG(SSPVEL(K)/SSPVEL(K-l)  )/GRADNT 
IF  (GRADNT  .EQ.  0.0)  T  -  T+(DEPTH(K)-DEPTH(K-1) )/SSPVEL(K) 

5  CONTINUE 

GRADNT  «  (HMV2 -SSPVEL (FINAL) )/(ZRCVR  -  DEPTH(FINAL) ) 

IF  (GRADNT  .NE.  0.0)  T  «  T  +  LOG (HMV2 /SSPVEL ( FINAL )) /GRADNT 
IF  (GRADNT  .EQ.  0.0)  T  -  T  ♦  (ZRCVR-DEPTH(FINAL)  )/HMV2 
BMVEL  >  ABS( (ZRCVR-ZPING)/T) 

ENDIF 

IF  (INTE  .EQ.  1)  THEN 
TEMP  -  ZPING 
ZPING  >  ZRCVR 
ZRCVR  ■  TEMP 
ENDIF 
RETURN 
END 


HYDEPTH  SUBROUTINE  TO  CALCULATE  TRANSDUCER  DEPTHS 


SUBROUTINE  HYDEPTH  (ZCABLE,  ZBUOY,  NUMSSV,  DEPTH, 

♦  SSPVEL,  HMVEL,  BTIME,  PHONEZ) 

CC  SUBROUTINE  TO  CALCULATE  THE  TEPTH  OF  THE  HYDROPHONE  BELOW  THE 
CC  BUOY. 

REAL*4  ZCABLE,  ZBOOY,  DEPTH,  SSPVBL,HMVEL,BTZME,ZOLO, PHONEZ 
INTEGER* 4  NUMSSV 

DIMENSION  DEPTM(40),  SSPVSL(40) ,PHONZZ(4) ,  BTIME(4,4) , ZBUOY(4) 
bo  5  I  -  1,  4 
ZB  *  ZBOOY  (I) 

ZOLO  -  ZCABLE 
, DO  3  KOUNT  ■  1,  10 

CALL  HMVELOC(ZB,  ZOLO,  NUMSSV,  DEPTH,  SSPVEL,  HMVEL) 

ZmW  -  HMVEL  •  BTIME(Z,I)  *  ZB 

IF  (ABS(ZNEW-ZOLO)  .LT.  1.0)  GO  TO  4 
ZOLO  ZNEH 

3  CONTINUE 

4  PHONEZ(I)  -  ZNEW 

5  CONTINUE 
RETURN 
END 


•  RPHONE  SUBROUTINE  TO  CALCULATE  RANGES  BETWEEN  TRANSDUCERS* 

SUBROUTINE  RPHONE  (RANGE,  HTIME,  PHONEZ,  NUMSSV,  DEPTH,  SLPVEL, 
♦  HMVEL) 

CC 

CC  RANCE(l)  -  PHONE  1  TO  PHONE  2 

CC  RANCE(2)  -  PHONE  3  TO  PHONE  2 

CC  RANCC<3)  -  PHONE  1  TO  PHONE  4 

CC  RANCE(4)  ■  PHONE  3  TO  PHONE  4 

CC  RANCE(S)  >  PHONE  2  TO  PHONE  4 
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CC  HANGE(6)  «  PHONE  1  TO  PHONE  3 

REAL*4  RANGE,  HTIME,  PHONEZ,  DEPTH,  SSPVEL,HKVEL,T 
INTEGER  NUMSSV 

DIMENSION  RANGE(6),  HTIKE(4,4),  PHONEZ(4),  DEPTH(40), 

*  SSPVEL(40) 

DO  1  I  »  1,  4 

1  CONTINUE 

T  -  (HTINE(2,I)  -  HTIME(1,2)/3.0)*1.5 

CALL  HMVELOC( PHONEZ (1),  PHONEZ (2),  NUMSSV,  DEPTH,  SSPVEL, 

*  HMVEL) 

RAMGE(l)  >  T  •  HMVEL 

T  «  0.5*HTIME(2,3)  0. S*HTIMB(3,2) 

CALL  HMVELOC( PHONEZ (2),  PHONSZ(3),  NUMSSV,  DEPTH,  SSPVEL, 

*  HMVEL) 

RANGB(2)  •  T  •  HMVEL 

T  ■  O.S*HTIME(l,4)  -t-  O.S*HTIME(4,l) 

CALL  HMVELOC(PHONEZ(l),  PH<M(BZ(4),  NUMSSV,  DEPTH,  SSPVEL, 

*  HMVEL) 

RAHGE(3)  -  T  •  HMVEL 

T  ■  (HTIMB(3,4)  -  HTIMB(4,3)/3.0)*1.5 

CALL  HMVELCC(PHONEZ(3),  PHONEZ(4),  NUMSSV,  DEPTH,  SSPVEL, 

*  HMVEL) 

RANGE(4)  -  T  •  HMVEL 

T  >  0.7S*HTIMS(2,4)  +  0.25*HTIMB(4,2) 

CALL  HMVELOC(PH(M(EZ(2),  PH(MIEZ(4),  NtmSSV,  DEPTH,  SSPVEL, 
>  HMVEL) 

RAMGE(S)  -  T  *  HMVEL 
T  •  0.2S*HTIMB(1,3)  ♦  0.7S*HTIMS(3,1) 

CALL  HMVELOC(PHONEZ(l),  PHOMEZ(3),  NUMSSV,  DEPTH,  SSPVEL, 

*  HMVEL) 

RANGE(6)  •  T  •  HMVEL 
RETURN 

END 


TNOT  SUBROUTINE  TO  CALCULATE  TRANSDUCER  POSITIONS  * 

»•••*•••••••••••••••••••••••••••••••••••*••••••••*•*•••••••• 

SUBROUTINE  TWOO  (RANGE,  XP,  YP,  ZP,  BUOYY) 

I  RXAL*4  RANGE,  BUOYY 

REAL*4  XP,  YP,  ZP,  OX,  DY,  ANS,  XSQRl,  XSQRS, 

*  YSQRl,  YSQR3,  RSQRl,  RSQRS,  RHS,  DISCRIN,  THETA,  PHI 
DIMENSION  BUOYY(4),  XP(4),  YP(4),  ZP(4),  RANGE(6), 

ANS(2,2) 

XP(3)iiSQRT(RANCE{6)**2  -  (SP(1)-ZP(3) )**2) 

YP(3)-0.0 

DX-XP(I)-XP(3) 

DY-YP(1)-YP(3) 

RSQRl-RANGE ( 1 ) **2- ( ZP ( 1 ) -ZP(2 ) ) **2 

RSQR3-RANGE(2)**2>(ZP(2)>SP(3))**2 

XSQR1-XP(1)**2 

XSQR3-XP(3)**2 

YSQR1-YP(1)**2 

YSQR3«YP(3)**2 

RHS  •  (XSQRl  *  YSQRl  •  XSQR3  >  YSQR3  -fRSQRl  -  RSQRl) *0.5 
IP  (DY  .EQ.  0.0)  THEN 
THETA  -  RHS/DX 
A  •  1 

B  -  -2  *  YP(1) 

C  «  (THETA)**2  ♦  XSQRl  ♦  YSQRl  -  (2  *  THETA  •  XP(1))  -  RSQRl 
DISCRIM  -  B*»2  -  (4  *  A  •  C) 
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ANS(1,2)  »  {-B  +  SQRT{ABS(DISCRIM) ) )/{2*A) 

ANS(2,2)  -  (-B/A  -  ANS{1,2)) 

ANS(1,1)  ■  THETA 
AMS(2,1)  *  THETA 
ELSE 

THETA  «  RHS/DY 
PHI  -  DX/DY 
A  -  1.0  +  PHI**2 

B  •  2*(YP(1)*PHI-THETA*PHI-XP(1)) 

C  -  XSQRl  +  (THETA)**2  +  YSQRl  -  RSQRl  -  2.0  *  THETA  *  YP(1) 
DISCRIM  -  P  *2  -  (4  *  A  *  C) 

ANS(1,1)  «  ,-B  ■^SQRT(DISCRIM))/(2*A) 

ANS(2,1)  -  '-B/A  -  ANS(1,1)) 

ANS(1,2)  «  THETA  -  PHI*ANS<1,1) 

ANS(2,2)  -  THETA  -  PHI*ANS(2,1) 

ENDIF 

IP  (ABS(ANS(1,2)  -  BUOYY(2)), .LT.  1000)  THEM 
XP(2)  «  ANS<1,1) 

YP(2)  -  ANS(1,2) 

VT  CV 

XP(2)  -  ANS(2,1) 

YP(2)  «  ANS(2,2) 

ENOIP 

OX-XP(2)-XP(3) 

DY-YP(2)-YP(3) 

RSQR1-RANGE(5)**2-(ZP(2)-ZP(4))**2 

RSQR3-RANGE(4)**2-(ZP(3)*ZP(4))**2 

XSQR1»XP(2)**2 

XSQR3«XP(3)**2 

YSQR1-YP<2)**2 

YSQR3»YP(3)**2 

RHS  ■  (XSQRl  ♦  YSQRl  -  XSQR3  -  YSQR3  -fRSQR3  -  RSQRl) *0.5 
IP  (DY  .SQ.  0.0)  THEN 
THETA  -  RHS/OX 
A  •  1 

B  •  -2  •  YP(2) 

C  •  (THETA)**2  ♦  XSQRl  ♦  YSQRl  -  (2  *  THETA  *  XP(2))  -  RSQRl 
DISCRIM  «  B**2  -  (4  *  A  •  C) 

ANS(1,2)  -  (-B  ♦  SQRT(DISCRIN))/(2*A) 

ANS(2,2)  -  (-B/A  -  ANS(1,2)) 

ANS(1,1)  -  THETA 
ANS(2,1)  -  THETA 
ELSE 

THETA  •  RHS/DY 
PHI  ■  DX/DY 
A  -  1.0  ♦  PHI**2 

B  ■  2*(YP(2)*PHI-THETA*PH1-XP(2)) 

C  ■  XSQRl  ♦  (THETA)**2  ♦  YSQRl  -  RSQRl  -  2.0  *  THETA  •  YP(2) 
DISCRIM  ■  B**2  -  (4  •  A  •  C) 

ANS(1,1)  •  (>B  •^SQRT(ABS(DZSCRIN)))/(2*A) 

ANS(2,1)  P  (-B/A  •  ANS(l,i;) 

ANS(1,2)  -  THETA  •  PHI*ANS(1,1) 

ANS(2,2)  •  THETA  -  PHI*ANS(2,1) 

ENDIP 

IP  (ABS(ANS(1,2)  •  BUOYY(4))  .LT.  ioOO)  THEN 
XP(4)  •  ANS(1,1) 

YP(4)  p  ANS(1.2) 

ELSE 
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XP(4)  -  ANS(2,1) 
yP(4)  -  ANS(2,2) 
ENDIF 
RETURN 
END 


•  TGTRANGE  SUBROUTINE  TO  CALCULATE  RANGE  TO  TARGETS 


SUBROUTINE  TGTRANGE  (TTIME.TGTZ,  ZPHONE,NUMSSV, DEPTH, SSPVEL, 
HMVEL,  TRANGE,  J) 

REAL*4  TGTZfZPHONE, DEPTH, SSPVEL,HMVEL,TRANGE,TTIME 
INTEGER*4  J,  NUMSSV 

DIMENSION  TGTZ ( 2 ) , ZPHONE ( 4 ) , DEPTH( 40) , SSPVEL ( 40 ) , TRANGE (2,4), 

♦  TTIME(2,4) 

DO  10  I  »  1,  4 

CALL  HMVELOC(TGTZ(J) ,  ZPHONE(I),  NUMSSV,  DEPTH,  SSPVEL,  HMVEL) 
TRANGE(J,I)  -  TTIME(J,I)  *  HMVEL 
10  CONTINUE 
RETURN 
END 


•  THREED  SUBROUTINE  TO  CALCULATE  TARGET.  POSITIONS 


SUBROUTINE  THREED  (J,TGTX,TGTy, TGTZ, PHONEX, PHONEY, PHONEZ,RM) 
REAL'*4  TGTX ,  TGTY ,  TGTZ ,  PHONEX ,  PHONEY  ,  PHONEZ ,  PRPX ,  PRPY ,  PRPZ  , 

*  XTO,YTO,ZTO,Al,A2,A3,A4,AS,A6,A7,A8,A9,RO,DR,RM 
ZNTEGER*4  J,  IN 

DIMENSION  TGTX(2),TGTY(2),TGTZ(2),PHONEX(4),PHONEY(4), 
PHONEZ(4),RO(4),RM(2,4) 

IN  -  0 
S  CONTINUE 
XTO-TCTX(J) 

YTOi^TY(J) 

ZTO«TCTZ(J) 

Al-0.0 

A2-0.0 

A3-0.0 

A4>0.0 

r.s-o.p 

*.6«0.6 

A7-0.0 

ABaO.O 

A9-0.0 

DO  10  I«  1,  4 

RO(I)  -  SQRT((XTO-PHONEX(I))**2+(YTO-PHONEY(I))**2+ 

♦  (ZTO-PHONEZ(I))**2) 

DR  -  RM<J,Z)  >  RO(Z) 

PRPX  >  (XTO  r  PHONEX(I))  /  RO(I) 

PRPY  -  (YTO  -  PHONEY(Z))  /  RO(I) 

PRPZ  >  (ZTO  -  PHONEZ(Z))  /  RO(Z)  . 


A1 

m 

A1 

♦ 

PRPX*PRPX 

A2 

m 

A2 

♦ 

PRPX*PRPY 

A3 

m 

A3 

■f 

PRPX •PRPZ 

A4 

m 

A4 

♦ 

PRPY*PRPY 

AS 

m 

AS 

♦ 

PRPY*PRPZ 

A6 

m 

A6 

♦ 

PRPZ*PRPZ 

A7 

m 

A7 

♦ 

DR* PRPX 

A8 

m 

A8 

♦ 

DR* PRPY 
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A9  «  A9  +  DR*PRPZ 
10  CONTINUE 

DBTAO  »  A1*(A4*A6-A5*A5)-A2«(A2*A6-A5*A3)+A3*(A2*A5-A4*A3) 
DETAl  ■  A7*{A4*A6-A5*A5)-A8*(A2*A6-A5*A3)+A9*(A2*A5-A4*A3) 
DETA2  ■  A1*(A8*A6-A9*A5)-A2*(A7*A6-A9*A3)+A3*(A7*A5-A8*A3) 
DETA3  ■  A1*(A4*A9-A5*A8)-A2*<A2*A9-A5*A7)+A3*(A2*A8-A4*A7) 
TGTX(J)  ■  XTO  +  DETAl /DETAO 

TGTY(J)  »  YTO  +  DETA2 /DETAO 

TGTZ(J)  -  ZTO  +  ’^1TA3/DETA0 

IP  (( (DETA1**2+DETA2**2+DETA3**2) /DETAO)  .LT.  1.0  )  THEN 
GO  TO  900 
ENDIF 

IN  -  IN  +  1 
IP  (IN  .GT.  25)  THEN 
GO  TO  900 
ENOIP 
GOTO  5 
900  RETURN 
END 


*  OUTPUT  SUBROUTINE  TO  OUTPUT  DATA  RESULTS 


SUBROUTINE  OUTPUT  (MEAN1,VAR1,MEAN2,VAR2,MEANR,VARR, 
+  BTN,STN,TTN,NCOUNT) 

REAL*4  MEAN1,VAR1,MEAN2,VAR2,MEANR,VARR,BTN,STN,TTN 
INTEGER*4  NCOUNT 

OPEN  (UNIT  •  8,  PILE  -  'OOTPOT3.DAT') 


WRITE(8,*) 

WRITE (8, 1000) 

WRITE (8, 1001) 

WRITE (8, 1002) 

WR1TE(8,*)  ' 

WRITE(8,1003)  NCOUNT,  BTN,  STN,TTN 


MEANl,  VARl 
MEAN2,  VAR2 
MEANR,  VARR 

NUMBER  OP  ITERATIONS,  BUOY 


MEAN 


STAND  DEV 


1000  FORMAT 

1001  FORMAT 

1002  FORMAT 

1003  FORMAT 
RETURN 
END 


(' 

(' 

( 


TRUE  ERROR  FOR  TARGET  1  » 
TRUE  ERROR  FOR  TARGET  2  - 
RELATIVE  ERROR  - 


(  3X,I6,15X,  F8.4,  F8.4,  F8.4) 


'»F8.3, 

',F8.3, 

',F8.3, 


TRANS 

F8.4) 

F8.4) 

F8.4) 


TARGET 


URAND  SUBROUTINE  TO  GENERATE  RANDOM  NUMBERS 


* 


SUBROUTINE  URAND (A, B, SEED) 
REAL*8  SEED,  PI,  Rl,  R2 
PI-3.14159265358979 
Rl*(SEED+Pl)**504.D-2 
Rl-Rl-DINT(Rl) 

SEEO-Rl 

R2-(SEED+PI ) **504. D-2 
R2-R2-DINT(R2) 

SSBD-R2 

R-SQRT ( -2  *LOO ( Rl ) ) 

T-2*PI*R2 

A-R*SIN(T) 

B-R*COS(T) 

RETURN 

END 
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APPENDIX  B 


The  following  is  the  output  data  from  the  simulations 
conducted  for  each  of  the  four  cases  discussed  in  Chapter  IV. 
Case  number  1: 


BUOYLOC.DAT 


0.0, 

0.0, 

5.0, 

0.0, 

0.0, 

1000.0 

3000.0, 

5196.0, 

5.0, 

3000.0, 

5196.0, 

1000.0 

6000.0, 

0.0, 

5.0, 

6000.0,. 

0.0, 

1000.0 

9000.0, 

5196.0, 

5.0, 

9000.0, 

5196.0, 

1000.0 

TARGET.DAT 


2036.0,  2036.0,  400.0, 
1965.0,  1965.0,  400.0, 

OUTPUT.DAT 

MEAN 

STAND  ; 

TRUE  ERROR  FOR  TARGET  1  - 

.002 

.0000 

TRUE  ERROR  FOR  TARGET  2  » 

.001 

.0000 

RELATIVE  ERROR 

.001 

.0000 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0000  .0000  .0000 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  «  19.829  .0122 

TRUE  ERROR  FOR  TARGET  2  -  19.695  .0122 

RELATIVE  ERROR  -  14.982  .0084 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 

1000  .0010  .0010  .0010 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  -  39.791  .0247 

TRUE  ERROR  FOR  TARGET  2  -  39.564  .0248 

RELATIVE  ERROR  -  30.101  .0169 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 

1000  .0020  .0020  .0020 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  -  60.054  .0380 

TRUE  ERROR  FOR  TARGET  2  -  59.780  .0383 

RELATIVE  ERROR  -  45.545  .0260 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0030  .0030  .0030 
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TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN  STAND  DEV 
*  80.887  .0530 

=  80.616  .0537 

=  61.654  .0362 

BUOY  TRANS  TARGET 
.0040  .0040  .0040 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

=  102.065 
=  102.172 
=  78.586 


STAND  DEV 
.0688 
.0718 
.0486 


BUOY 

.0050 


TRANS 

.0050 


TARGET 

.0050 


‘PRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
REX  ATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

=  123.875 
=  123.837 
=  98.655 


STAND  DEV 
.0910 
.0939 
.0726 


BUOY 

.0060 


TRANS 

.0060 


TARGET 

.0060 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

145.096 

146.658 

118.392 


STAND  DEV 
.1044 
.1133 
.0931 


BUOY 

.0070 


TRANS 

.0070 


TARGET 

.0070 


MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  «  164.285  .1207 

TRUE  ERROR  FOR  TARGET  2  -  165.526  .1290 

RELATIVE  ERROR  -  138.519  .1075 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0080  .0080  .0080 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  -  183.286  .1335 

TRUE  ERROR  FOR  TARGET  2  «  179.448  .1264 

RELATIVE  ERROR  »  154.021  .1197 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0090  .0090  .0090 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  -  198.663  .1406 

TRUE  ERROR  FOR  TARGET  2  ■  199.931  .1478 

RELATIVE  ERROR  -  172.593  .1288 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0100  .0100  .07.00 
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Case  number  2 


BUOYLOC.DAT 

0.0,  0.0,  5.0, 

3000.0,  5196.0,  5.0, 

6C./0.0,  0.0,  5.0, 

9000.0,  5196.0,  5.0, 

TARGET.DAT 

1000.0,  2000.0,  400.0, 
3000.0,  4236.0,  400.0, 

OUTPUT.DAT 

TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
REIATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TOUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


0.0,  0.0,  1000.0 
3000.0,  5196.0,  1000.0 

6000.0,  0.0,  1000.0 
9000.0,  5196.0,  1000.0 


MEAN  STAND  DEV 

-  .001  .0000 

»  .002  .0000 

a:  .001  .0000 

BUOY  TRANS  TARGET 

.0000  .  .0000  .0000 

MEAN  STAND  DEV 

*»  23.314  .0147 

-  12.332  .0058 

»  19.779  .0111 

BUOY  TRANS  TARGET 

.0010  .0010  .0010 

MEAN  STAND  DEV 

»  46.881  .0300 

»  24.694  .0116 

-  39.783  .0226 

BUOY  TRANS  TARGET 

.0020  .0020  .0020 

MEAN  STAND  DEV 

-  71.010  .0467 

•  37.106  .0175 

«  60.300  .0355 

BUOY  TRANS  TARGET 

.0330  .0030  .0030 

MEAN  STAND  DEV 

-  96.323  .0673 

»  49.592  .0236 

•  81.926  .0521 

BUOY  TRANS  TARGET 

.0040  .0040  .0040 
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TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN  STAND  DEV 
=  122.966  .0913 

=  62.178  .0299 

=  104.788  .0721 

BUOY  TRANS  TARGET 

.0050  .0050  .0050 

MEAN  .  STAND  DEV 
=  148.598  .1136 

=  74.547  .0364 

=  127.339  .0906 

BUOY  TRANS  TARGET 

.0060  .0060  .0060 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

»  174.871 
=  87.047 
=  150.839 


STAND  DEV 
.1362 
.0428 
.1110 


BUOY  TRANS  TARGET 
.0070  .0070  .0070 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

194.483 

99.075 

169.387 


STAND  DEV 
.1489 
.0498 
.1218 


BUOY  TRANS  TARGET 
.0080  .0080  .0080 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

209.383 

111.316 

183.449 


STAND  DEV 
.1497 
.0572 
.1207 


BUOY  TRANS  TARGET 
.0090  .1.090  .0090 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN 

232.755 

123.187 

206.640 


STAND  DEV 
.1665 
.0622 
.1412 


BUOY  TRANS  TARGET 

.0100  .0100  .0100 


32 


Case  nvunber  3 


BUOYLOC. 

DAT 

0.0, 

0.0, 

5.0 

3000.0, 

5196.0, 

5.0 

6000.0, 

0.0, 

5.0 

9000.0, 

5196.0, 

5.0 

TARGET.DAT 

2036.0,  2036.0,  400.0, 
1965.0,  1965.0,  400.0, 


OUTPUT.DAT 

TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


0.0,  0.0,  1000.0, 

3000.0,  5196.0,  900.0, 

6000.0,  0.0,  800.0, 
9000.0,  5196.0,  700.0, 


MEAN  STAND  DEV 

-  .002  .0000 

«  .001  .0000 

.001  .0000 

BUOY  TRANS  TARGET 

.0000  .0000  .0000 

MEAN  STAND  DEV 

»  21.890  .0139 

»  21.668  .0138 

«  16.342  .0095 

BUOY  TRANS  TARGET 

.0010  .0010  .0010 

MEAN  STAND  DEV 

44.028  .0283 

•  43.639  .0283 

•  32.941  .0193 

BUOY  TRANS  TARGET 

.0020  .0020  .0020 

MEAN  STAND  DEV 

«  66.771  .0445 

«  66.268  .0440 

•  50.214  .0303 

BUOY  TRANS  TARGET 

.0030  .0030  .0030 

MEAN  STAND  DEV 

«  90.263  .0629 

«  89.303  .0620 

»  68.667  .0436 

BUOY  TRANS  TARGET 

.0040  .0040  .0040 
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TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  I'RROR 

NUMBER  C'F  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE,  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN  STAND  DEV 

=  114.024  .0834 

*  114.166  .0862 

«  89.502  .0640 

BUOY  TRANS  TARGET 

.0050  .0050  .0050 

MEAN  STAND  DEV 

*  137.825  .1066 

=  136.027  .1048 

*  111.802  .0864 

BUOY  TRANS  TARGET 

.0060  .0060  .0060 

MEAN  STAND  DEV 

=*  164.585  .1340 

«  159.147  .1289 

=  136.492  .1275 

BUOY  TRANS  TARGET 

.0070  .0070  .0070 

MEAN  STAND  LEV 

=«  178.254  .1318 

=  174.554  .1306 

*  146.033  .1085 

BUOY  TRANS  TARGET 

.0080  .0080  .0080 

MEAN  STAND  DEV 

»  201.941  .1519 

«  198.101  .1477 

=  178.767  .1563 

BUOY  TRANS  TARGET 

.0090  .0090  .0090 

MEAN  STAND  DEV 

»  217.521  il554 

212.118  .1509 

»  186.845  .1463 

BUQV  TRANS  TARGET 

.0x00  .0100  .0100 


Case  number  4:, 


BUOYLOC.bAT 


0.0, 

,  0.0, 

5.0, 

0.0, 

0.0, 

1000.0 

3000.'), 

5196.0, 

5.0, 

3000.0, 

5196.0, 

900.0 

oOUO.O. 

0.0, 

5.0, 

6000.0, 

0.0, 

800.0 

9000.0, 

5196.0, 

5.0, 

9000.0, 

5196.0, 

700.0 

TARGET.DAT 

1000.0,  2000.0,  400.0, 

3000.0,  4236.0,  400.0, 

OUTPUT.DAT 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  =  .001  .0000 

TRUE  ERROR  FOR  TARGET  2  =  .001  .0000 

RELATIVE  ERROR  =  .001  .0000 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 

1000  .0000  .0000  .0000 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  =  24.415  .0156 

TRUE  ERROR  FOR  TARGET  2  *  13.655  .0068 

RELATIVE  ERROR  =  20.218  .0114 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 

1000  .0010  .0010  .0010 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  »  49.162  .0319 

TRUE  ERROR  FOR  TARGET  2  -  27.372  .0136 

RFLATIVE  ERROR  *  40.721  .0233 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 

1000  .0020  .0020  .0020 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  »  74.694  .0504 

TRUE  ERROR  FOR  TARGET  2  =  41.198  .0208 

RELATIVE  ERROR  »  61.921  .0372 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0030  .0030  .0030 

MEAN  STAND  DEV 

TRUE  ERROR  FOR  TARGET  1  *  102.194  .0764 

TRUE  ERROR  FOR  TARGET  2  •  55.192  .0283 

RELATIVE  ERROR  -  84.955  .0585 

NUMBER  OF  ITERATIONS,  BUOY  TRANS  TARGET 
1000  .0040  .0040  .0040 
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TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


TRUE  ERROR  FOR  TARGET  1 
TRUE  ERROR  FOR  TARGET  2 
RELATIVE  ERROR 

NUMBER  OF  ITERATIONS, 
1000 


MEAN  STAND  DEV 

=  130.319  .1024 

=  69.366  .0364 

=  109.312  .0842 

BUOY  TRANS  TARGET 

.0050  .0050  .0050 

MEAN  STAND  DEV 

-  154.690  .1222 

=  83.145  .0447 

=  130.410  .0952 

BUOY  TRANS  TARGET 

.0060  .0060  .0060 

MEAN  STAND  DEV 

=  183.675  .1501 

=  97.235  .0547 

=  156.063  .1203 

BUOY  TRANS  TARGET 

.0070  .0070  .0070 

MEAN  STAND  DEV 

=  202.344  .1559 

=*  112.773  .0695 

=*  173.353  .1238 

BUOY  TRANS  TARGET 

.0080  .0080  .0080 

MEAN  STAND  DEV 

»  221.332  .1593 

«  126.445  .0794 

-  194.134  .1413 

BUOY  TRANS  TARGET 

.0090  .0090  .0090 

MEAN  STAND  DEV 

-  239.836  .1740 

»  140.192  .0861 

-  210.705  .1459 

BUOY  TRANS  TARGET 

^0100  .0100  .0100 
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